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ABSTRACT 

To investigate further the comparison between AGN and black hole X-ray binaries, we 
have studied the main X-ray variability properties of the Seyfert 1 Galaxy Mrk 335. 
We put particular emphasis on the X-ray time lags, which is a potentially important 
diagnostic of physical models. From a 100 ksec observation by XMM-Newton we show 
that the power spectrum of this source is well fitted by a bending power law model, 
and the bend time-scale Tb is precisely at the value predicted by the Tb vs H^ line- 
width relation of McHardy et al. Variations in different energy bands show time-scale 
dependent time lags, where higher energy bands lag lower ones. The lag, r, varies as 
a function of the Fourier frequency, /, of the variability component in the light curves 
as r oc f^ 1 at low frequencies, but there is a sharp cut-off in the lags at a frequency 
close to the bend frequency in the power spectrum. Similar behaviour is seen in black 
hole X-ray binary systems. The length of the time lags increases continuously with 
energy separation, in an almost log-linear relation. We show that the lag spectra can 
be produced by fluctuations propagating through the accretion flow as long as the 
energy spectrum of the X-ray emitting region hardens towards the centre. 

Key words: Galaxies: active 



1 INTRODUCTION 

Understanding the similarities between Active Galactic Nu- 
clei (AGN), which are powered by accretion onto super- 
massive black holes, and the much smaller Galactic black 
hole X-ray binary systems (BHXRB) is currently one of 
the major research topics in high energy astrophysics. Al- 
though some properties might be expected to scale simply 
with mass, the fact that the accretion disc is much cooler in 
AGN may affect properties which depend on the spectrum 
of the nuclear ambient photon field, e.g. any X-ray emission 
produced by Comptonisation. 

Over recent years considerable attention has been de- 
voted to the comparison of the X-ray variability properties 
of AGN and BHXRB and, in particular, to comparison of 
their X-ray power spectrum. It has been shown that most 
X-ray bright AGN vary in a manner similar to that of the 
persistent BHXRB Cyg X-l in the high/soft state (where the 
medium energy X-ray spectrum is soft and the flux is high) , 
rather than in the low/hard state (where the m edium energy 
X-ray spectrum is hard and the flux is low) (e.g. Uttlev et al.l 
I2002J ; iMcHardv et"afll2004 120051 ; lUttlev fc M c Hardvll2005lh 
In the soft state, the power density spectrum (PDS) of 
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Cyg X-l has a power law shape of slope -1 at low frequencies 
(i.e. P(f) oc /~ Q with a — 1), bending to a steeper slope 
(a Js 2) above a bend frequency /&, and the PDS of most well 
observ e d AGN are similar. W e have shown IMcHardv et al.l 
l|2007t ); [K"ording et ail (|2007l ') that the time-scale associated 
with the bend, Tb ~ l/fb, scales as black hole mass divided 
by accretion rate (in Eddington units) from AGN down to 
BHXRB. As higher quality AGN X-ray data has become 
available, primarily from XMM-JVewton, it is now possible 
to investigate other X-ray timing measurements which may, 
ultimately, prove to have greater diagnostic power than sim- 
ple bend time-scales, i.e. X-ray time lags and coherence. 

Time lags are normally observed between the fluctua- 
tions in different energy bands, where the hard X-ray bands 
tend to lag softer ones and the coherence is the degree to 
which the hard and soft bands are correlated. Where the 
data are good enough, the lags and coherence can be mea- 
sured as a function of Fourier frequency. A particularly im- 
portant aim, to which observations of lags and coherence can 
contribute, is to understand the origin of the X-ray varia- 
tions in both BHXRB and AGN. An interesting model for 



the or igin of the fluctuations was proposed by iLvubarskiil 
(1997). In this model, characteristic time-scales of variabil- 
ity, produced by variations in accretion rate, are associated 
with each radius in the accretion disc, with longer time- 
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scales being associated with larger radii, and the variations 
propagate inward towards the X-ray emitting region. This 
model provides a natural explanation of the observed lin- 
ear relationship be tween rms variability and flux in both 
BHXRB and AGN (|Uttlev fc M c Hardvll2001j ). 

The propagating fluctu ation model has been ex tended 
bv lKotovet"aiT(|200ll ) and I Arevalo fc Uttlevl [|2006ri to in- 
clude an extended X-ray emitting region where the X-ray 
spectrum hardens towards the black hole. As the fluctua- 
tions, which might have been produced further out, travel 
through this emitting region, they modulated first the emis- 
sion with a softer spectrum and later the emission with 
harder spectra, thereby naturally giving rise to hard lags. 
The model has been used to simultaneously explain both 
the lags and the X-ray PDS in the hi gh accretion rate Nar - 
row Line Seyfert 1 (NLSf) Ark 564 (JMcHardv et al.ll2007T l 
but it has not yet been widely applied to other AGN. In 
this paper we analyse a long XMM-Newton observation of 
the X-ray bright NLS1 Mrk 335 and show how the prop- 
agating fluctuation model can simultaneously explain the 
lags and the PSD. Using our new measurement of Tb, we 
also investigate how well Mrk 335 agrees with the relation- 
ship between Tb and the width of the Hp emission line of 
JMcHardv et alj|2007l l. 

The paper is organised as follows: we briefly discuss the 
data reduction in section [2] calculate the power spectrum 
and fit it with a bending power law model in Sec. [3] We 
study the relation between different energy bands in Sec. [3] 
for the coherence and Sec.[5]for the time delays. In Sec.[6]we 
summarize the main variability properties of Mrk 335 found 
in this work. 



2 THE DATA 

Mrk 335 was observed by XMM-JVewton for ~ 130 ks in 
2006 January 3-5 during revolution 1112. Spectral fitting 
and spectral variability analyse s of these data have bee n 
published bv lO'Neill et all J2007h and lLarsson et al.1 l |2007h . 

We used dat a from the EPIC PN detector 
(jStriider et al.l |2001[ ). The PN camera was operated in 
Small Window mode, using thin filter. The data were 
processed using XMM-SAS v 6.5.0. Source photons were 
extracted from a circular region of 5' in radius and a 
background region of equal area was chosen on the same 
chip. Source and background events were selected by quality 
flag=0 and patterns=0-4, i.e. only single and double events 
were used. The background level was generally low and 
stable, except at the beginning and at the end of the 
observation. We discarded the first 10 ks and last 12.4 ks to 
obtain 110.6 ks long light curves free of background flares. 

For the power spectrum and cross spectrum analyses, 
we constructed light curves in the 0.2-0.6, 0.6-1, 1-3 and 
3-10 keV bands, using 54 sec long time bins, with average 
count rates of 9.4, 4.5, 4.2 and 0.8 counts/sec, respectively. 
The bin size was chosen to obtain 2 11 points in each light 
curve, which enables us to use the fast Fourier transform 
(FFT) in the power spectrum fitting procedure. 




bending power law model fit 
to the 0.2-10 keV PDS 
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Figure 1. 0.2-10 keV PDS of Mrk 335 plotted in blue dots with 
error bars, with the best-fitting bending power law model plotted 
as a solid black line. The power due to Poisson noise, which starts 
to dominate the variability at 10 — 3 Hz, has been subtracted. 



3 POWER SPECTRUM 

All the measured AGN power spectra have an approxi- 
mately P oc 1// shape, breaking to a steeper slope at high 
frequencies (e.g lUttlev et al.1 12002|; IMarkowitz et al.l 120031 ; 
iMcHardv et al.l 120051 . |2004| ; iMarkowitz et all 120071 '). similar 
to the PDS of Cyg X-l in the high/soft state. In only one 
AGN, Ark 564, an additional, low frequency break has been 
significantly de tected, below which t he PDS levels off to a 
constant value (|McHardv et al.ll2007l ). This shape resembles 
Cyg X-l in either the low/hard or very high states, which 
have roughly a doubly-broken power law shape. In either 
case, single or double break, the high frequency end of the 
PDS has approximately a bending power-law shape with a 
low frequency slope of -1. For the PDS of Mrk 335 we will 
consider a low frequency slope of this value to test whether 
a high frequency bend is required by the data within the 
observed frequency band. 

We constructed the PDS of Mrk 335 l ight curves 



through the discrete Fourier Transform (DFT) (|Press et al.l 
[1993), 



P(fi. 



2 At 



Re 2 x (fi) + lm 2 x (fi)\, 



(1) 



where Kex(fi) and Inix(/;) are the real and imaginary parts 
of the discrete Fourier transform of the time series x(t), x 
is the average count rate, At is the sampling time interval 
and N is the number of points in the light curve. With this 
normalisation, the power spectrum is given in terms of a 2 /x 2 
with units of 1/Hz so the integral of P(f) over / equals the 
normalised variance. We binned the periodogram points in 
frequency bi ns, adopting the logari t hmic binning method 
described bv lPapadakis fc Lawrence! |l993j). A minimum of 
2 points were included in each frequency bin and the bin 
width was set at 1.3 x/i, where /; is the smallest Fourier 
frequency in the corresponding bin. 

We fitted a bending power law model defined as 



P(f) 



Af~ 



(2) 



to the PDS of the 0.2-10 keV light curve. The low frequency 
slope oll = 1 was fixed and the high frequency slope an, 
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Figure 2. Confidence contours of the bending power law fit to 
the 0.2-10 keV PDS. The 68%, 90% and 99% levels are shown by 
the solid, dashed and dotted lines, respectively. The best fitting 
values are marked by the black dot. 



the bend frequency /;, and the normalisation A were left 
as free parameter s. We used the Mon te Carlo fitting proce- 
dure PSRESP of lUttlev et all lJ2002h . to take into account 
possible distorting effects of red-noise leak and aliasing and 
produce robust error bars and goodness of fit criteria. 

The parameter space was searched over a grid of mul- 
tiplicative factor 1.5 in /(,, from 10~ 5 to 10 -2 ' 5 Hz, and an 
additive factor 0.2 in an, from 1 to 4. A finer grid, with 
steps of 1.1 in /(, and 0.1 in an was used to search over the 
best fitting region and the resulting parameter values were 
essentially the same. 

The best-fitting bending power law model, together 
with the data PDS, unfolded by observational biases, is 
shown in Fig. [l] This model produces a good fit with ac- 
ceptance probability of 87%. The high frequency slope has 
a value of an — 3.8 with a 90% confidence lower limit of 2.5 
and a 99% limit of 2.2. This slope is therefore significantly 
different to the low frequency slope of 1 and the bend in the 
PDS is clearly detected. The best-fitting bend frequency has 
a value of 1.7^12 x 10 _4 Hz, where the uncertainties corre- 
spond to ler errors. The lower 90% confidence limit on the 
bend frequency, however, falls below 10 -5 Hz, which is the 
lowest frequency probed. Figure [2] shows the confidence con- 
tours of this fit over the entire parameter space explored. 
As shown below, the bend in the PDS appears to shift to 
higher frequencies for higher energy bands. The 99% con- 
fidence lower limit on the bend frequency falls inside the 
available parameter space for the 1-3 and 3-10 keV bands. 
Figure [3] shows the confidence contours of the fit to the 1-3 
keV band PDS, where the bend frequency is well bounded. 

Fitting the 0.2-10 keV band PDS, allowing the low fre- 
quency slope ai to vary produces essentially the same val- 
ues: an = 3.8^", fb = 1.6_ 12 x 10~ 4 Hz , and a best fitting 
value of cv; = 1.1, but no useful constraints can be placed 
on CV;. 



3.1 Bend frequency vs Ha relation 

iMcHardv et all J20o3) have shown for a sample of AGN with 
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Figure 3. Confidence contours of the bending power law fit to 
the 1-3 keV PDS. The 68%, 90% and 99% levels are shown by 
the solid, dashed and dotted lines, respectively. The best fitting 
values are marked by the black dot. In this energy band and also 
in 3-10 keV, the lower 90% and 99% confidence limits on the 
break frequency are contained in the parameter space explored, 
i.e. at frequencies higher than the lowest bin frequency which is 
determined by the length of the observation. 
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Table 1. PSRESP fits to the PDS in different energy bands using 
a bending power law model with fixed low frequency slope a; = 1. 
The errors correspond to 66% confidence contours in the ajj vs 
/(, plane, asterisks denote unbounded contours. 



fb measured as above, that fb is well correlated to the width 
of the H/3 line in the optical spectrum, following the relation 



logTL = 4.2 log #,3,fwhm - 14.43, 



(3) 



where TJ, is the break time-scale in units of da y s and 
-H/3,fwhm is in units of km/s. iPeterson et al.l (|2004f l 
used data from two observational campaigns to measure 
H^fwhm of Mrk 335 obtaining values of 1629T145 and 
1375T357 km/s. According to Eq. [3] these line-widths pre- 
dict a break time-scale of 0.078 < Tb < 0.16 days in the 
first case and 0.016 < Tb < 0.15 days in the second. The 
measured bend frequency corresponds to a time-scale of 
0.068l J 2 3 days, consistent with both line width measure- 
ments, within the uncertainties. Figure [3] shows the position 
on Mrk 335 on the l ine-width vs bend time-scale plot of 
IMcHardv et al.l (120071 ) . for both line- width measurements. 



3.2 Energy dependence of the PDS 

The PDS of Mrk 335 shows more high frequency power 
at higher energies, as is normally observed in AGN (e.g. 



McHardv et alj2004 iMarkowitz et al.ll2007l ; IMcHardv et all 
2007h . Fitting the 0.2-0.6, 0.6-1, 1-3 and 3-10 keV PDS 
with PSRESP, keeping the low slope fixed to 1, results in 
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Figure 4. H,3 line width vs break time-s cale: the solid black line 
represents the best-fitting relation from iMcHardv et al.l (2007) 
who used the 9 AGN plotted as red circles and the low luminosity 
AGN NGC4395 plotted as an open crossed circle. The 2 AGN 
(green filled squares) whose upper limit on the break time-scale 
is not constrained are plotted, but were not included in the fit. 
The black crosses represent the position of Mrk 335 on this plot, 
using the best-fitting break time-sca le of the 0.2-10 keV P DS and 
the two line width measurements of lPeterson et al.l (2004). 



higher bend frequencies for higher energy bands as shown in 
Table [T] although they are partly within the error from each 
other. The errors on these parameters are large, however, 
because they take into account the fact that the observation 
is only one realisation of a large set of possible observed PDS 
with the same underlying power spectrum. Considering that 
the observations in the different energy bands were made 
simultaneously, the uncertainty in the relative PDS shape 
is much smaller. To illustrate the energy dependence more 
clearly, we calculated the ratios of the harder band PDSs to 
the 0.2-0.6 keV PDS, which is shown in Fig. [5] If the energy 
bands were completely coherent, the only source of scatter 
would be the Poisson noise, which adds different amounts 
of power to the PDS points of each light curve, while the 
stochastic deviations of the coherent variations cancel out. 
The variations in different bands are not completely coher- 
ent, as shown below, so we considered the reduced coher- 
ence when calculating the error bars on the PDS ratio. We 
assumed that 10%, 20% and 40% of the fluctuations in the 
0.6-1, 1-3 and 3-10 keV light curves is incoherent with the 
0.2-0.6 keV fluctuations and added the uncertainty in this 
incoherent part of the power to the uncertainty in the power 
of the Poisson noise. These levels of incoherent fluctuations 
were chosen to match the average coherence functions be- 
tween the corresponding energy bands, at least at low fre- 
quencies where the coherence can be well measured, as de- 
scribed in Sec. [4] below. With these error bars, the PDS 
ratios are not consistent with a simple change in high fre- 
quency slope to account for the additional high-frequency 
power, as in that case the ratios would increase as a power 
law instead of leveling-off and dropping. Equal power law 
slopes and different bend frequencies produce PDS ratios of 
a stepped shape, that rise rapidly at frequencies above the 
lower bend frequency and level-off above the higher bend fre- 
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Figure 5. Additional high-frequency power at high energies: 
Power spectra in the 0.6-1 (black), 1-3 (red) and 3-10 keV 
(blue) bands divided by the PDS of the 0.2-0.6 keV band. Higher 
energy bands show increasingly more variability power above 
~ 2 X 10~ 5 Hz. 



quency. To produce a drop at even higher frequencies would 
require a steeper high-frequency slope for the hard band. 

For visualization purposes, we over-plotted in Fig.[5]the 
ratios between bending power law model PDSs, with low fre- 
quency slope oti = 1 and different bend frequencies and a^. 
For the model lines shown, we used the best fitting /t, from 
TableQ]for the corresponding energy bands and fixed ah — 3 
for the soft band. We allowed aj, of the hard bands to vary in 
order to fit the amplitude of the PDS ratios of the data. The 
best fitting values of a^ were 3.5 for 0.6-1 keV, 3.6 for 1-3 
keV and 4.8 for 3-10 keV. The shape of the PDS ratio indi- 
cates that higher energy bands have both a higher bend fre- 
quency and a steeper high frequency slope. This behaviour 
could indicate the presence of a band-limited high frequency 
variability component with a hard energy spectrum, adding 
on to the soft band PDS. This energy dependence, however, 
can also be produced by extended emitting regions filtering 
out the high frequency variability power. If higher energies 
are emitted by a more compact region, e.g. if the emitted 
spectrum hardens towards the centre, their PDSs will have 
higher bend frequencies than softer bands, while the high 
fr equency slope can be eq ual or steeper (see figures 2 and 6 
in lArevalo fc Uttlevll2006r i . 



4 COHERENCE 

We will now investigate the correlation between variations in 
different energy bands. The coherence -y 2 , i.e. the degree of 
linear correlation between two light curves, can be computed 
as a function of Fourier frequency using the cross spectrum 
C(/) = S* (/)#(/), where S(f) and H(f) are the Fourier 
transforms of the soft and hard band light curves s(t) and 
h(t), respectively. 



i 2 Ui) = 



(Re c (/0> 2 + (Im c (/ i )> 2 



(4) 



^S(/0l 2 X|H(/0l 2 > 

where Rec(/0 and Imc(/0 are the real and imaginary parts 
of the cross spectrum C(f) and angle brackets represent av- 
eraging over independent measurements. In our case, the in- 
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dependent measurements are taken from consecutive Fourier 
frequencies in a frequency bin. Coherence values range from 
1, f or totally coherent , t o for totally incoherent variations. 
See IVaughan fc Nowakl (|l997l ) for a detailed description of 
the coherence function interpretation and error calculation. 

Coherence functions were calculated between 0.2-0.6 
keV and each of the three harder bands (0.6-1, 1-3 and 
3-10 keV) and the results are shown in Fig. [6] A minimum 
of 8 points was included in each frequency bin and the bin 
separation was fixed at 1.5 x/j. 

As the Poisson noise and other distorting effects in the 
observed light curves can artificially decrease the coherence, 
we applied the Poi s son n oise correction factor described in 
IVaughan fc Nowakl (1 19971 ) and estimated any additional de- 
viations using simulated light curves. For each pair of energy 
bands, we generated 1000 p airs of coherent red-nois e light 
curves using the method of iTimmer fc Konid (|l995), with 
the underlying PDS parameters obtained from the fits to the 
corresponding real light curves. The average count rates of 
the simulated light curves were made to match the measured 
values and Poisson deviates were added to each simulated 
light curve point. The coherence was then calculated for each 
simulated pair using the same sampling and Poisson noise 
correction method as used for the real data. The mean of the 
distribution of simulated coherence measurements for each 
frequency bin is plotted as a red solid line in Fig. [6] and 
the top and bottom 95% extremes of this distribution are 
plotted as dashed blue lines in the same plot. 

This analysis shows that the variations in the 0.2- 
0.6 keV and 0.6-1 keV energy bands are highly coherent 
(7 ~ 0.95) up to frequencies of at least 2 x 10~ 4 Hz. Above 
this frequency Poisson noise makes the coherence estimates 
unreliable so it is not possible to detect a real drop in co- 
herence. At larger separations between the energy bands the 
coherence drops significantly, even at the lowest frequency 
probed, 4 x 10~ 5 Hz, to ~ 0.87 between 0.2-0.6 and 1-3 keV 
and to ~ 0.6 between 0.2-0.6 and 3-10 keV light curves. 



5 TIME LAGS 

Time delays between variations in different energy bands 
can be measured using the phase lag spectrum </>(/). This 
spectrum gives the difference in phase between the Fourier 
transform components of the soft and hard light curves, as a 
function of Fourier frequency /, and it can be c onverted into 
a time -lag spectrum as r(/) = </>(/)/(2-7r/). See lNowak et al.l 
(1999) for a full description of the lag spectrum measure- 
ment. In what follows we will refer to the time lag spectrum 
simply as lag spectrum. 

In all AGN where the lag spectrum has 
been measured, the hard bands lag softer ones 
(Papadakis, Nandra & Kazanas 2001; Vaughan et al. 2003; 



McHardv et all 12004 lArevalo et al.ll2006l ; llVIarkowitz et all 
20071 ) and the size of the lag increases with decreasing 
Fourier frequency. In the BHXRB Cyg X-l in the high/soft 
state, the lag spectrum has a power law dependence on 
frequency with a slope of ~ 0.7 and a cut-off at high 
frequencies, close to the break in the PDS. In the low/hard 
state, the lag spectrum shows approximately the same 
power law dependence but with additional structure in the 
form of steps. AGN data are normally not good enough 



to detect such structure and only simple power laws with 
fixed slope are fitted. A change in slope in the lag spectrum 
of an AGN has only been detected significantly in Ark 564 
|Arevalo et al.l 120060 but in that case the 'step' was too 
deep to resemble the lag spectrum of Cyg X-l in the 
low/hard state and rather resembles a Very High state. 



5.1 Lag spectrum 

The lag spectrum of Mrk 335 between 0.2-0.6 and 3-10 keV 
energy bands is shown in Fig. [7] where positive values indi- 
cate hard band lagging. A minimum of 2 points per bin was 
used and the phase-lag measurements were averaged over 
frequency bins of width 2x/;, before being converted into 
time-lags. Notice that the lags are shown only for frequen- 
cies below 3x 10 _4 Hz, where the coherence is still acceptably 
high. 

As seen in other AGN and XRB, the hard band in 
Mrk 335 lags the softer band and the lag size decreases 
with increasing frequency with an approximately power law 
dependence. The slope, however, appears steeper than the 
t(/) oc / _ model normally fitted to AGN lag spectra. 

We used simulated light curves to determine whether 
the steepness of the slope is produced by observational bi- 
ases. Light curve pairs were generated as in Sec. [4] intro- 
ducing time lags between the simulated soft and hard light 
curves of the form r(/) = 0.082/~ 1 , where the amplitude of 
the lag spectrum was obtained from a fit to the real data. 
The simulated light curves were sampled following the pat- 
tern of the real data and lag spectra were calculated for each 
pair using the same method and binning. The mean of the 
distribution of 1000 simulated pairs is shown by the red dot- 
ted line in Fig. [7] and the top and bottom 95% extremes of 
the distribution are plotted in blue dashed lines. The real 
lag spectrum, plotted in black markers joined by a solid line 
in the same figure, drops more steeply than the simulations, 
indicating that the drop in high-frequency lags is not due to 
observational biases on an intrinsically 1// lag spectrum. We 
repeated the same experiment using the best-fitting power 
law model of the real data, r(f) = 8.8 x 10~ 6 /" 184 , as the 
underlying lag model for the simulations. In this case, the 
real data falls within the 95% extremes of the distribution of 
simulated lag spectra, but this model over predicts the scat- 
ter, putting almost the entire 95% lower limit on negative lag 
values. Fitting the lag spectrum with a range of slopes be- 
tween -1 and -1.84 and repeating this experiment produced 
simulated lag distributions that are either too tight and do 
not include the high-frequency lag points and/or too wide 
and over predict the scatter in the measured lag spectrum 
at low frequencies. Therefore, the lag spectrum is not consis- 
tent with a single powe r law and the bend do wn or cut-off at 
around 10 -4 Hz is real. lArevalo et al.l (12006) found a similar 
change in slope in the lag spectrum of Ark 564, in that case 
however, the bend was found a decade in frequency below 
the bend frequency of the PDS, while in Mrk 335 the bend 
in the PDS and lag spectrum appear at the same frequency, 
within the uncertainties. 



5.2 Energy dependence of the time-lags 

The length of the time lags increases with the separation 
between the energies of the bands considered. We calculated 
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Figure 6. Coherence as a function of Fourier frequency between 0.2-0.4 and 0.4—0.6 keV energy bands (left), 0.2-0.6 and 1-3 keV 
bands (centre) and 0.2—0.6 and 3—10 keV bands (right), plotted in solid black lines. The dot-dashed red lines represent the median of 
a distribution of simulated coherence functions, constructed from intrinsically coherent light curves affected only by Poisson noise. For 
the simulations we used the count rate and underlying PDS shape of the corresponding light curves. The dashed blue lines represent the 
95% upper and lower limits of this distribution, marking the scatter expected in the case of perfectly coherent light curves, distorted by 
Poisson noise and the sampling pattern in the same way as the real data. Evidently, the measured coherence is significantly lower than 
1 between the 0.2-0.4 and 1-3 keV bands and between the 0.2-0.4 and 3—10 keV bands, even at the lowest frequencies sampled. 
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Figure 7. Time lag as a function of Fourier frequency / between 
the 0.2—0.4 keV and 3—10 keV energy bands over the frequency 
range where the coherence is well measured, the data points are 
shown in black markers joined by the solid line. The dotted and 
dashed lines represent the mean and top and bottom 95% limits of 
a distribution of simulated lag spectra with an underlying lag = 
0.082/ -1 lag spectrum (see text). The drop of the lags at high 
frequencies is greater than expected for this underlying model, 
showing that the real lag spectrum has cither a slope steeper 
than -1 or it cuts-off at high frequencies. 



lag spectra of the 0.2-0.4 keV to the 0.4-0.6, 0.6-1, 1-2, 2- 
5 and 5-10 keV bands, finding that the higher the energy 
band, the more it lags the 0.2-0.4 keV band. The fractional 
lags, r/time-scale, are approximately 0.72 % for 0.4-0.6 keV, 
1.5% for 0.6-1 keV, 4.2% for 1-2 keV, 6.2% for 2-5 keV and 
12% for 5-10 keV. Each lag spectrum resembles the one 
shown in Fig. [7] having a cut-off at high frequencies. These 
fractional lag values only represent the low-frequency part 
of the spectra. 

In Fig. [8] we plotted the lowest- frequency lag for each 
pair as a function of the ratio between the average energy 
of the bands compared. The lags increase in an almost log- 
linear manner with increasing energy ratio. The best fitting 




Figure 8. Time lag at frequency 1.3 X 10~ 5 Hz as a function of 
energy. The logarithm of the ratio between the average energies 
of the each pair of energy band is plotted in the x axis. The solid 
line represents a linear fit to the data plotted. 



log-linear model, shown by the solid line in the plot has a 
form 



r= (8317 ±1138) log 



-Eh a 



BL 



(1585 ± 580) 



(5) 



For a perfect log-linear relation, the intercept should equal 
zero, i.e. lags between identical energy bands should be 0. 
The intercept of the fitted relation is significantly different 
from this value so the behaviour of the time lags in Mrk 33 5 
is not entirely consistent with Cyg X-l (jNowak et al.lll999l ). 
Notice however, that the log-linear relation in Cyg X-l was 
calculated for energies above 2 keV and in Mrk 335 we used a 
baseline energy band of 0.2-0.4 keV. It is possible that other 
spectral components in Mrk 335 below 2 keV, in particular 
the soft excess might be suppressing the intermediate-energy 
lags. 



5.3 Interpretation of the lag spectrum 

We applied the propagating fluctuation paradigm of 
iLvubarskiil (|1997l ) to attempt to reproduce the observed lag 
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spectra. In this scenario, the fluctuations are produced by 
the accretion flow over a large range in radius, and travel to- 
wards the centre to modulate the X-ray emitting region. We 
assume that the fluctuations are produced and propagated 
on the viscous time-scale of a thick accretion flow and that 
the emitted energy spectrum hardens towards the centre, 
therefore producing hard lags. 

The travel time tt of a fluctuation of frequency /, from 
its radius of origin Tf to any radius r is given by 

/ dr 



"M = / zh 



v{r) 3(H/R) 2 a 



('/ 



3/2 3/2x 



(6) 



where the second equality assumes a propagation speed 
v(r) — (H/R) 2 ar~ ' 2 in units of c and r in units of 
R g = GM/c 2 . This speed is taken fro m the standard disc 
prescription (jShakura fc Sunvaevlll973r ). where a is the vis- 
cosity parameter and H/R the disc thickness to radius ratio. 
The time lags arise from different radial emissivity pro- 
files e(r) characterising the different energy bands. The 
weighted travel time, f (/) of a fluctuation on frequency / 
folded through a given emissivity profile is 



r(/) = 



JJ 7 tt(r)e(r)2Trrdr 



(7) 



fj e(r)2ivrdr 

We will replace the normalisation factor J f e(r)rdr by E in 
what follows and cancel the 2-7T factors. Finally, the time lag 
between two energy bands with emissivity profiles t a (r) and 
et (r) is the difference between the corresponding weighted 
travel times, 



r(f) = fa(f)-Mf) = 



ttlr 



gq(r) 

E a 



tt(r) 



e b (r) 



E b 



rdr(8) 



replacing tt{r) by the expression in Eq. [6] and cancelling 
terms gives: 



"(/) = 



3(H/R) 2 a 



7 „5/2 f eajr) _ e b {r) 

E a Eb 



dr 



(9) 



Notice that the length of the time lag is not simply pro- 
portional to 1/(H/R) 2 a because the integration limit r/ is 
also a function of these disc parameters. If the disc thick- 
ness or a decrease, the propagation speed decreases making 
the lag longer, but the fluctuations on a given frequency are 
produced at a smaller radius r/ , so they travel a shorter dis- 
tance and so t(/) remains at approximately the same level. 

Equation [5] gives the lag spectrum between energy 
bands characterised by emissivity profiles e a and et- From 
standard disc theory, the total energy dissipated per unit 
area follows a e ex r -3 relation. As a simple example we 
used similar power law emissivity profiles but with higher 
exponents for higher energy bands, to produce a spectrum 
that hardens towards t he ce ntre. 

iKellv fc Bechtoldl (|2007h measured a mass for Mrk 335 
of 6.6^32 x 1O 6 M0, we used this mass estimate to calculate 
Rg/c = 33sec to compare model predictions to the real data. 
In this model, both the bend in the PDS and the cut-off 
in the lag spectrum are related to the characteristic time- 
scale of the innermost region of the accretion flow. This nor- 
mally requires thick disc parameters to produce the power 
observed at the highest frequencies. We used (H/R)a = 0.1 
to put the viscous time-scale, t v — r/v v i sc (r), at a radius of 
6R g (i.e. the last stable orbit of a non-rotating black hole) 
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Figure 9. Time lag as a function of Fourier frequency between the 
0.2-0.4 keV energy band and consecutively higher energy bands. 
The length of the lag increases with the separation of the energy 
bands. The dashed lines represent model fits as explained in Sec. 
5.3 



at 10 4 s for a 6.6 x 10 6 Mq black hole, which is appropriate 
for modelling Mrk 335. 

The measured lag spectra between the 0.2-0.4 keV band 
and all harder bands are shown by the markers and solid 
lines in Fig. \§\ where the higher energy hard bands produce 
longer lags. The quality of the data does not allow a detailed 
fit to each lag spectrum but we attempted a simple imple- 
mentation of the model described above to make a broad 
comparison. The dashed lines in Fig. [5] represent the model 
fits, where we assumed parameters (H/R)a — 0.1, inner- 
most radius of the disc ri n = 6R g and an emissivity profile 
for the 0.2-0.4 keV band of e oc r fi with f3 . 2-0.4 = -3.0. 
The other energy bands were also assumed to have power 
law emissivity profiles and the exponents were allowed to 
vary in order to fit the data. The standard viscous propaga- 
tion speed resulted too high and produced lags consistently 
shorter than the observed values. We reduced the propaga- 
tion speed by a factor of two for the fits, which reproduced 
the observed average lags. The models shown in Fig. [9] have 
exponents of -3.08, -3.15, -3.33, -3.6 and -3.8 for the 0.4-0.6, 
0.6-1, 1-2, 2-5 and 5-10 keV bands, respectively. 

The model does reproduce the power law shape of each 
lag spectrum and predicts a cut-off at the right frequency. 
The range of lag spectra amplitudes can be reproduced with 
moderate differences in the emissivity profiles of the energy 
bands. These energy-dependent emissivity profiles can be 
translated into a radial dependence of the emitted spec- 
trum. The exponents stated above produce spectral slope 
changes from Ar = T(r-2) — I\ri) = 0.35 log^/fi) when 
we compare 0.2-0.4 to 0.4-0.6 keV energy bands up to 
Ar = 0.61og(r 2 /ri) when comparing 0.2-0.4 to 5-10 keV 
energy bands, where T2 > n, flux(i?) oc E~ v and E is the 
average energy of the band. This difference in spectral hard- 
ening implies that, for the model to work exactly, either the 
radial emissivity profiles are not exact power laws or that 
the emitted spectrum at each radius is not a power law. 

As Mrk 335 is accreting at a high rate, an accretion 
disc with constant sc ale height H, at least a t the c entre, can 
be more appropriate. ISvensson fc Zdziarskil |l_994J) solve the 
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structure equations of a thin disc plus corona system, and 
show that the accretion flow represented by the corona is 
radiation pressure dominated at accretion rates above ap- 
proximately 18% the Eddington limit. As this is probably 
the case in Mrk 335, if we identify the corona with the ac- 
cretion flow producing and propagating the variability, we 
should consider a scale-height radial dependence of H/R = 
r _1 m(l — y/6/r)/2. Assuming that the same viscous veloc- 
ity prescription applies in this case, then the propagation ve- 
locity will contain an additional factor of r~ 2 which reflects 



the change in (H/R) with radius, i.e. 



(H/R) 2 ar 



-5/2 

where (H/R) 2 , is a fixed value. The resulting travel time of 
the fluctuations changes to tt(r) — 7(H ? R yi ( r J — T ), 
which can then be used in Eq. 8 to calculate the lag spectra. 
The shape of the lag spectra are approximately the same as 
in the case of constant (H/R) 2 , shown in Figure [9] As long 
as the factor (H/R) 2 a is the same at the innermost radius, 
the constant H/R and constant H models produce the bend 
at the same frequency. The constant H model however pro- 
duces even shorter lags, by approximately a factor of 4, so 
for this model to produce the observed lag values the propa- 
gation speed should be reduced by a factor of ~ 8 compared 
to the viscous velocity. 



6 DISCUSSION AND CONCLUSIONS 

We have studied the standard variability properties of 
Mrk 335 using a 114 ks observation by XMM-Newton. The 
main results are summarized below. 

• The Power spectrum is well described by a bending 
power law model, with low frequency slope fixed at ai — 1, 
best-fitting bend frequency of ft, — 1.7l?;| 5 x 10~ 4 Hz and 
high frequency slope of a h = 3.8t** . 

• The measured bend frequency coincides precisely with 



the value predi c ted by the ft, vs H,3 line-width relation of 
iMcHardv et al.l 1120071). using the H^ line-width measure- 
ments of IPeterson et al.l (|2004j), which demonstrates the va- 
lidity of this relation. 

• The power spectrum is energy-dependent with more 
high-frequency power at higher energies. The difference in 
the PDSs of different energy bands indicates a shift in the 
PDS bend to higher frequencies for higher energy bands. 

• The degree of coherence between variations in different 
energy bands decreases with the energy separation of the 
bands considered. The coherence is significantly lower than 
1, between the 0.2-0.6 and the 1-3 keV and 0.2-0.6 and 3-10 
keV energy bands, even at the lowest frequencies sampled 
(4xl0~ 5 Hz). 

• Time lags are found between variations in different en- 
ergy bands, where harder X-rays lag softer bands. The lag 
spectrum as a function of Fourier frequency can be fit with 
a single power law of slope ~ —1.8, much steeper than the 
lag spectra observed in other AGN. Alternatively, the lag 
spectrum of Mrk 335 is consistent with a power law of slope 
~ — 1 dropping to a steeper slope at around 10 _4 Hz, where 
the bend in the PDS is also found. The latter shape resem- 
bles the high frequency end of lag spectra of Cyg X-l in both 
low/hard and high/soft states. 

• The size of the time lags increases with energy sepa- 
ration following an almost log-linear relation. In all cases 



tested, using the energy bands 0.2-0.4, 0.4-1, 1-2, 2-5 and 
5-10 keV, any band lags all softer ones, at least at low fre- 
q uencies where the e rr or bar s are small. 

lArevalo fc Uttlevl l|2006l ) have shown tha t the 
propagating-fluctuation paradigm of iLvubarskiil (119971 ) 
reproduces the shape of the PDS and the rms-flux relation. 
This variability model is successful in producing variability 
over a very broad range of time-scales and it naturally 
predicts a cut-off at high frequencies. Complemented with 
an extended emitting region whose spectrum hardens 
towards the centr e, this model can explain the time lags 
l|Kotov et al.ll200ir ) and is compatible (but does not ensure) 
high coherence in the fluctuations of different energy bands. 
The extent of the emitting region in the model also acts as 
a low-pass filter reducing the high frequency power in the 
more extended energy bands, i.e. the softer bands, shifting 
the bend frequency to lower values. 

The energy dependence of the PDS, together with the 
shape of the lag spectra of Mrk 335 are consistent with the 
predictions of this model, if the propagation speed is reduced 
by a factor of 2-8, compared to the viscous velocity, depend- 
ing on the disc structure. The lag spectra can be reproduced 
by a simple implementation of the model, with an emitted 
energy spectrum that hardens towards the centre. The en- 
ergy dependence of the length of the time lags requires a 
maximum hardening of the spectrum with decreasing ra- 
dius, between any radii r\ and r^, of Ar = 0.61og(r2/ri). 
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